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Gauge dependence of the fermion quasiparticle poles in hot gauge theories 
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The gauge dependence of the complex fermion quasiparticle poles corresponding to soft collective 
excitations is studied in hot gauge theories at one-loop order and next-to-leading order in the high- 
temperature expansion, with a view towards going beyond the leading order hard thermal loops and 
resummations thereof. We find that for collective excitations of momenta k ~ eT the dispersion 
relations are gauge independent, but the corresponding damping rates are gauge dependent. For 
k <C eT and in fc — > limit, both the dispersion relations and the damping rates are found to be 
gauge dependent. The gauge dependence of the position of the complex quasiparticle poles signals 
the need for resummation. Possible cancellation of the leading gauge dependence at two-loop order 
in the case of QED is briefly discussed. 

PACS numbers: ll.lO.Wx, 11.15.-q, 12.38.Cy 



I. INTRODUCTION 

The demonstration of gauge independence of phys- 
ical quantities is of fundamental importance in gauge 
theories at finite temperature. Interestingly, the main 
breakthrough in perturbative gauge theories at high tem- 
perature was in fact motivated by the quest for gauge- 
independent damping rates of plasma excitations. It was 
realized [1] that in hot gauge theories the usual connec- 
tion between the order of the loop expansion and powers 
of the gauge coupling constant is lost, hence contributions 
of leading order in the gauge coupling constant arise from 
every order in the loop expansion. Resummation to all 
orders of these leading, gauge-independent hard thermal 
loops (HTLs) [1-6] is necessary for a consistent calcu- 
lation. The HTL resummation program [1, 2] leads to 
an effective theory that systematically includes contribu- 
tions of different momentum scales [7] . 

It was proved in Ref. [8] that the singularity struc- 
ture (i.e., the position of poles and branch singularities) 
of gauge boson and fermion propagators are gauge inde- 
pendent when all contributions of a given order of a sys- 
tematic expansion scheme are accounted for. While the 
position of quasiparticle poles in the leading order HTL 
approximation are completely gauge independent [9, 10], 
nevertheless gauge dependence will enter at subleading 
order and gauge-independent extensions beyond the lead- 
ing order HTL results are not yet available. More re- 
cently, several authors [11-13] have shown that the trun- 
cated on-shell two-particle-irreducible (2PI) effective ac- 
tion has a controlled gauge dependence, with the explicit 
gauge-dependent terms always appearing at higher order. 
It would be interesting to study possible cancellation of 
the leading gauge dependence in the singularity struc- 
ture of gauge boson and fermion propagators beyond the 
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leading order HTL approximation. 

The propagation of fermions in hot and dense matter 
is of interest in a wide variety of physically relevant sit- 
uations. In ultrarelativistic heavy ion collisions and the 
formation and evolution of the quark-gluon plasma, lep- 
ton pairs play a very important role as clean probes of 
the early, hot stage of the new state of matter [14]. The 
propagation of quarks during the nonequilibrium stages 
of the electroweak phase transition is conjectured to be 
an essential ingredient for baryogenesis at the electroweak 
scale in (non)supersymmetric extensions of the standard 
model [15]. In stellar astrophysics, electrons and neutri- 
nos play a major role in the evolution of dense stars such 
as white dwarfs, neutron stars and supernovae [16]. 

The goal of this article is to investigate the gauge de- 
pendence of the complex fermion quasiparticle poles cor- 
responding to soft collective excitations at one-loop or- 
der and next-to-leading order in the hight-temperature 
expansion. There have been investigations of how the 
leading HTL dispersion relations of fermions (and the 
gauge independence thereof) are affected by retaining 
nonleading powers of temperature [17, 18]. While it is 
known that quasiparticles at finite temperature will in 
general acquire thermal widths due to coUisional broad- 
ening rendering the quasiparticle poles complex, the pre- 
vious authors considered mainly the real parts of the 
complex quasiparticle poles (dispersion relations) with- 
out discussing the corresponding imaginary parts (damp- 
ing rates). In this article we study the gauge depen- 
dence of both the real and imaginary parts of the com- 
plex quasiparticle poles on equal footing. Furthermore, 
we focus on the soft fermion collective excitations of mo- 
menta k < eT, where e <C 1 is the gauge coupling con- 
stant and T the temperature. On the one hand, it allows 
us to fill gap in the literature where either a numerical 
analysis with k ^ T was carried out [17] or the limiting 
cases oi eT ^ k ^ T and fc = were considered [18]. 
On the other hand, k < eT is the relevant momentum re- 
gion at the heart of the HTL resummation program [1, 2]. 
Progress made in the soft momentum region will certainly 
shed light on the issue of going beyond the HTL resum- 



mat ion. 

The rest of this article is organized as fohows. In Sec. II 
we calculate the real and imaginary parts of the one- loop 
fermion self-energy in general covariant gauges and up to 
next-to-leading order in the high-temperature expansion. 
In Sec. Ill we study the gauge dependence of the complex 
fermion quasiparticle poles corresponding to soft collec- 
tive excitations of momenta k < eT. Possible cancella- 
tion of the leading gauge dependence at two-loop order in 
the case of QED is briefly discussed. Sec. IV presents our 
conclusions. In the Appendix we calculate the vacuum 
contribution to the fermion self-energy that is neglected 
in the main text. 



II. ONE-LOOP FERMION SELF-ENERGY IN 
COVARIANT GAUGES 

We will carry out our perturbative calculations in 
the imaginary-time (Matsubara) formalism (ITF) of fi- 
nite temperature field theory [9, 10]. The continuation 
to imaginary time that describes the theory at finite 
temperature T is obtained by replacing it -^ t with 
< T < /3 = 1/T. Contrary to the usual ITF, however, 
we will not work in Euclidean spacetime but keep the 
metric tensor and the Dirac gamma matrices the same 
as in Minkowski spacetime. 

In ITF it proves convenient to work in the spectral rep- 
resentation of the propagators. The fermion propagator 
is given by 



s(zi/„,k) 



p/(A:o,k) 
dko-r : — , Vn 



(2n+l)7rT, (1) 



where pf is the free Dirac fermion spectral function (zero 
chemical potential) [9] 

P/(fco,k)=^sgn(fco)'5(i^'), if=(fco,k), (2) 

with sgn(x) being the sign function. The gauge boson 
propagator is given by 



#"(*^„,p) 



dpo 



P^"(P0,P) 



ujn — 2mrT, (3) 



where p'^'^ is the free gauge boson spectral function in 
general covariant gauges [19] 



P^^bcP) = sgn(po) 



(i-^)P'P''l^ 



d_ 

dpi 



6{P^ 



(4) 

with P = (pojP)- The covariant gauge parameter ^ is 
defined in such a way that ^ = 1 is the Feynman gauge. 
Since at one-loop order the fermion self-energy has 
the same structure (up to gauge group factors) in both 
Abelian and non-Abelian gauge theories, for notational 
simplicity wc will consider the Abelian case in what fol- 
lows. The non-Abelian case can be obtained through the 
replacement e^ -^ g^Cp, where Cp = {N"^ — l)/2N is 



the Casimir invariant of the fundamental representation 
in SU(A^) gauge theories. 

In QED the one-loop fermion self-energy in ITF is 
given by 



E(w„,k) = e^T^ / — -g7^s(w„ + Wm,q)7'' 



xd^^(iw„,p). 



(5) 



where q = k -I- p. Upon substituting the fermion and 
gauge boson propagators into (5), the sum over the 
bosonic Matsubara frequency can be done easily [9] lead- 
ing to 

f (fp f°° f°° 

i;(z;/„,k) = eM —33 / dpn / rf(7opA"^(Po, p) 

^l^Pf{<la,Ci)Y : , (6) 

(Jo - PO - IVn 

where nB,F{x) ~ l/{e^^ =F 1) are the Bose and Fermi 
distribution functions, respectively. 

After the analytic continuation w„ -^ lo + 0^ to arbi- 
trary frequency u, the imaginary part of the (retarded) 
self-energy can be readily found to be given by 



ImE(a;, k) 



e TT 



d^p 



1,^1'j J-00 J-oc 

^l''Pfiqa,<lh''[nBiPo) + npigo)] 
xS{uj-qo+po). (7) 

The real part of the self-energy is obtained from the imag- 
inary one through the dispersive representation 



ReS(a;,k)==PV / — 

J-00 T^ 



dko ImS(A:o,k) 



ko ~ LU 



(8) 



where PV denotes the principal value. 

For massless fermion, rotational invariance and chiral 
symmetry entail that the fermion self-energy in equilib- 
rium can be parameterized by [9, 10] 

S(a;, k) = S(0) {u, k) 7O + S(i) (c^, fc) 7 • k, (9) 

where k = \k\, Vi = k/fc, and S*^°) and S^-'^' are scalar 
functions 

E(")(c.,fc) = itr[S(^,k)7"], 

E(i)(a>,fc) - -itrp(a>,k)7.k]. (10) 

Since we are interested in the self-energy calculated in 
general covariant gauges, it is convenient to decompose 
the former into gauge-independent and -dependent con- 
tributions. We write 



E(a;, k) — Sfg('-^i k) -|- S^(a;, k) 



(11) 



where Spc is the gauge-independent part calculated in 
the Feynman gauge (FG) and S^ is the remaining gauge- 
dependent part. 



A. The imaginary part 

We first study the gauge-independent contribution 
IniEpG- Evaluating the Dirac traces and performing the 
trivial integrals over pq and qq in (7), we obtain 



.(0) 



ImE^"(^(co,fc) 



e TT 



3 — {[1 + "s(p) - npiq)] 



lmJ:^pl{uj,k) = -e^TT -— : 



(27r)3 2p 

xS{p + q-uj) + [nB{p) +nF{q)] 
xS{p- q + uj)} + (lli ^ -uj), 
(Pp k • q 



(27r)3 ^p-{[^ + nB{p)-np{q)] 

xS{p + q-uj) + [uBip] +np{q)] 
x5{p-q + uj)}- {lj —> -Lj). (12) 

In the above expressions the support of the energy- 
conserving delta functions S{p + q — uj) and S{p + q + Lj) is 
to > k and u> < k, respectively, corresponding to the usual 
two-particle cuts, while that oi6{p—q^uj) is uj'^ < k"^ cor- 
responding to the Landau damping cut which is purely a 
medium effect [1, 9]. 

The different contributions in (12) have a physical in- 
terpretation in terms of (off-shell) scattering processes 
taking place in the thermal medium. The terms propor- 
tional to 5{p + q — uj) arise from the processes in which a 
(timelike) fermion decays into a fermion and a gauge bo- 
son /* -^ /+7, and those proportional to 5(p—q+bj) orig- 
inate in the Landau damping process in which a (space- 
like) fermion scatters off a gauge boson in the medium 
/* -I- 7 ^- /. Here the off-shell fermion is denoted by a 
superscript "*" . 

In what follows we will neglect the vacuum contribu- 
tion, i.e., the terms which do not contain any thermal 
distribution functions, as we are interested mainly in 
the finite-temperature medium effects. Interested readers 
can find in the Appendix the calculation for the vacuum 
contribution. 

The angular integration over 77 — k-p in (12) can 
be performed analytically by using change of variables 
rj -^ z — q{ri) and the energy-conserving delta functions. 
The requirement that the energy-conserving delta func- 
tions must have a nonempty support restricts the range 
of the radial integration over p for fixed k and lu. The 
remaining integration over p for arbitrary k and uj is an 
involved numerical task [17] which, however, is not very 
useful for the purpose of studying next-to-leading order 
corrections. In order to compare with the leading order 
HTL results, we here focus on a high-temperature expan- 
sion for which k,uj <^T and keep terms in ImE to 0{T) 
in the Landau damping cut contribution but to 0(T°) in 
the two-particle cut contribution. 

A comment here is in order. One would presumably 
expect that it is sufficient to keep terms in ImE uniformly 
up to next-to-leading order in the high-temperature ex- 



pansion, namely, terms of order 0{T). Whereas this is 
correct for generic fc ~ w <C T, it is not correct in the 
limit k <^ Lo <^ T which is in turn relevant to the effec- 
tive fermion mass and the damping rate at rest that is 
obtained in the long- wavelength limit {k —^ 0). As will 
be seen below, the ©(T^) two-particle cut contribution in 
ImE will give rise to 0{hiT) term in ReE which will be- 
come ncxt-to-leading order in T in the limit k <^ lo <^ T. 
Separating the leading HTL contributions that can be 
calculated analytically, then in the remaining contribu- 
tions expanding the distribution functions in the high- 
temperature limit and cutting off the potentially diver- 
gent momentum integrations at T, we obtain after some 
algebra {k,Lu < T) 



.(0) 



lmT,^p^{uj,k) 



16k 

xd{k^ 



1 



2lu 



TT^T 



■In 



LJ 



Sir 



uj — k 

LO 



In 



OJ 



-1- 



.(1) 



lml?^^{uj,k) 



X In 



H 

2T 

IQk 



e{u;^ 



U! — k 

LJ 1 

k^ 2 
k 
2T 



-sgn(a)) 



k 

e{k^- 
^ k^ 



In 



LO 
LO + k 



LO 



-k 



(13) 



where 0{x) is the Hcaviside step function. Three im- 
portant features are gleaned from the above expressions: 
(i) The gauge-independent part ImEpc receives contri- 
butions both from above and below the light cone, cor- 
responding respectively to the two-particle and Landau 
damping cuts, (ii) The leading 0{T'^) terms in ImEpc 
are recognized as the HTL results, which arise solely from 
the Landau damping process [1, 5, 6, 9]. (iii) The sub- 
leading terms are suppressed by inverse powers of T and 
originate in the Landau damping as well as in the fermion 
decay processes. 

Next, we calculate the gauge-dependent contribution 
ImE J. The gauge-dependent part of the free gauge bo- 
son spectral function contains the derivative of the on- 
shell delta function d5{pQ — p'^)/dpQ, thus the integral 
over po can be done using integration by parts which 
in turn gives rise to derivatives of the Bose distribution 
function dnsip) / dp as well as of the energy-conserving 
delta functions d5{p — q — bo)/dLo, etc. Such structures 
are expected to be generic to the gauge-dependent higher 
loop contributions linear in (^ — 1). 

After some lengthy but straightforward algebra we ob- 
tain 



.(0) 



lm^'>{LO,k) 



duB (p) 



lTJiE^P{uJ,k) 



4&^(1 + P-^) 



[1 + nsip) - nFiq)]-^5{p + q - w) 



+ 



dp 
dnsip) 



S{p + q- uj) 



dp 



S{p-q + u) 



- (1 - p • q; 



[nsip) + nF{q)]-Q^S{p - q + lu) 



2 / d^P f (k-p)(p-q)-k-q 



4(27r) 



P 



[l + nsip) -npiq)] 



xS{p + q-uj) + [riBip) + nFiq)]S{p - q + uj) 



kp 



(1 + p-q) 



X ni + neip) - nF{q)]-^S(p + q - lu) - '^^^^' S{p + q - lu) 

+ (1 - p • q)([nB{p) + nF{q)]-^S(p - q + lu) + ^^^j^6{p - q + lu) 
— (lu^ —lu). 



(14) 



It is noted that the various contributions in (14), be- 
ing gauge dependent and hence in contrast to those in 
(12), wiU not have a physical interpretation in terms 
of the scattering processes taking place in the medium. 
The gauge-dependent contributions are expected to can- 
cel in physical quantities in a consistent calculation that 
generally requires resummation of perturbation theory. 
Clearly such a task is beyond the scope of this article. 

The angular integration in (14) can be done analyt- 
ically as in the gauge-independent contribution. Again 
neglecting the vacuum contribution, we find in the high- 
temperature limit the following rather compact expres- 
sions {k,uj <C T) 



:f' 



ImEr^(u>,fc) 



Idtt 

-M0(,2_ 



1 






.(1) 



lmSy'{LU,k) ~ (C-1) 



El 



k') 

LU 



-k 



lu 

fc2 



In 



-sgn(w) -0(0)2 



UJ — k 
(15) 



where we have combined part of the contributions from 
above and below the light cone by using the identity 
9{x) + 0{—x) = 1. The leading gauge-dependent terms in 
IniE^ are of order 0(T) and the subleading ones are sup- 
pressed by inverse powers of T as is the case for IniEpG- 



B. The real part 

From the above results for the imaginary part of the 
fermion self-energy the corresponding real part can be 
obtained using (8). In the high-temperature limit that 



we are interested in, the upper and lower limits of the 
integral over ko in (8) can be cut off at T and — T, re- 
spectively. The contributions to ReS from the ignored 
regions of integration are at most of order 0{T^), hence 
are negligible at the next-to-leading order under consid- 
eration. It is worthy noting that since ImS'-^-' [ImE'^-'^'] is 
an even (odd) function of tj therefore, as a result of (8), 
ReS(°^ [ReS'^^] is an odd (even) function of lo. 

In the high-temperature limit we obtain for the gauge- 
independent part {k,LU <^T) 



Re4°(^(a),fc) 



Re4G(w,fc) 



' 16fc 



167r2 

8k 

e^k 

'l67r2 



■In- 



2fc 



In 



-k 



;2-fc2| 



In- 



2k 

rj.2 



2k 



-In 



fc2|' 



UJ — k 
(16) 



where we have kept terms up to ©(InT) as remarked 
above. This is because the 0{T) term in Rei;(°) [ReSf^)] 
is proportional to k (fc2) in the limit k <^ uj <^T , there- 
fore the putative subleading ©(InT) term becomes next- 
to-leading order in T in this limit [note that ReS*^^^ (a>, k) 
vanishes identically in the limit fc ^ due to rotational 
symmetry] . 

Again the leading 0{T'^) terms in RcSfg are the HTL 
results [3, 5, 6, 9, 10]. The ©(InT) terms in (16) arise 
from the 0{T^) terms above the light cone in (13), thus 
originating in the region of hard loop momentum. Such 
In T contributions are not unique in the one- loop fermion 
self-energy. Indeed, similar InT behavior that originates 
also in the hard loop momentum region has been found in 
the gauge boson polarization [4, 20] as well as, in general. 



the n-point vertex function [21] at one-loop order in high 
temperature non-Abehan gauge theories. 

Similarly we find for the gauge-dependent part (fc, lo <C 
T) 



.(0) 



Re^'\uj,k) ~ -(C-1) 



m ■ 



167r2 |t^2_fc2|' 



p2t T2 

ReE^'^(a>,fc) ~ (^-1)— ^In. 



.(1), 



167r2 |cj2-fc2|- 



(17) 



The InT contributions in (17), like the gauge- 
independent ones in (16), arise from the region of hard 
loop momentum as well. Following the argument in 
Ref. [5] based on the ultraviolet divergences in the ab- 
sence of distribution functions, one might expect the 
leading term in ReE^ to be linear in T. Nevertheless, 
explicit calculation shows that the leading term actually 
goes like InT at high temperature. As we will see mo- 
mentarily, the absence of the 0{T) term in ReE^ have 
an important consequence in the gauge dependence of 
fermion dispersion relations. 

Furthermore, upon comparing (16) and (17), we find 
that the O(lnT) term in ReE^ has the same prefactor as 
that in the corresponding ReEpc: hence they can be com- 
bined together to yield a single gauge-dependent term 
proportional to ^. This, however, is not the case for the 
imaginary part of the self-energy, where only some of the 
terms in ImE^ and ImEpG can be combined together. 



III. GAUGE DEPENDENCE OF THE FERMION 
QUASIPARTICLE POLES 

Having calculated the one-loop fermion self-energy at 
next-to-leading order in T in general covariant gauges, we 
now proceed to study gauge dependence of the complex 
fermion quasiparticle poles. 

The full (retarded) inverse fermion propagator is given 
by the Dyson-Schwinger equation 



iS'^iij, k) = tJ7° - A:7 • k + E(tj, k). 



(18) 



where E(ti;,k) is the (retarded) fermion self-energy. 
Eq. (18) can be inverted to yield [9, 10] 



5(^,k) 



7° - 7 • k 7° + 7 • k 
A+(a),fc) ^ A_(w,fc) 



where 



with 



A±(a;, k) — uj^ k + E±(a;, fc), 



-v(0) 



E± (c^, fc) == E^^^ (t^, fc) ± E^^^ (cj, fc) 



(i)r 



(19) 



(20) 



(21) 



The analytic continuation of the fermion propagator to 
complex frequency features the following singularities: 



• Isolated poles. Isolated real poles of the fermion 
propagator correspond to stable quasiparticle exci- 
tations, whereas complex poles to unstable excita- 
tions (resonances) with finite widths. 

• Branch cuts. Branch cuts of the fermion propaga- 
tor correspond to multiparticle states. 

As we are interested in the collective excitation in the 
medium, we will consider isolated poles in the rest of 
this section. 

Write uj = E — i^ with E and 7 real. In the nar- 
row width approximation for which 7 <C -E, the equa- 
tions that determine the position of the complex poles 
are given by [22] 



E^k + KeY.±{E,k) = 0, 
7 + sgn(7)Z±(fc)ImE±(E,/c): 



0, 



(22) 



where Z±{k) in the second equation are the residues at 
the poles (wave function renormalizations) 



Z±{k) 



aReE±(w,fc) 



(23) 



u=E 



If the product of the residue at the pole and the imag- 
inary part of the self-energy on the quasiparticle mass 
shell Z±{k)I'mJ^±{E,k) is negative, there are two com- 
plex poles conjugate to each other in the physical sheet 
corresponding to a growing and a decaying solution, i.e., 
an instability. On the other hand, if the product is pos- 
itive there is no complex pole in the physical sheet, the 
pole has moved off into the unphysical (second or higher) 
sheet. In this case the spectral function features a Breit- 
Wigner shape resonance with a width given by the damp- 
ing rates 



7(fc) = Z±(fc)ImE±(£;,fc), 



(24) 



which in general determines an exponential falloff e"''''^'^^ * 
of the fermion propagator in real time. It is worth not- 
ing that numerical results reveal that at next-to-leading 
order the product Z±(fc)ImE±(_E, fc) is positive in the 
Feynman gauge. 

From the one-loop fermion self-energy calculated 
above, we recover at leading order in T [i.e., 0{T'^)] the 
celebrated HTL results [3, 5, 6, 9, 10]. There are two 
branches of stable fermion collective excitations (fermions 
and the so-called plasminos [9, 10]) with positive and neg- 
ative helicity to chirality ratios, respectively. Their dis- 
persion relations are manifestly gauge independent and 
given by (for positive energy solutions) 



E{k) 



a;+(fc) fermion, 
a;_(fc) plasmino. 



(25) 



A plot of the HTL dispersion relations uj±{k) can be 
found in the literature (see, e.g., Refs. [9, 10]). One of 
the main features is that the two branches of collective 
excitations develop a gap m = eT/vS, corresponding to 



an effective thermal mass that respects both gauge in- 
variance and chiral symmetry [3, 5, 6, 9, 10]. 

The gauge dependence of the quasiparticlc poles at 
next-to-leading order must be studied with care because 
the 0{T) terms in ReE(a;, k) becomes subleading in the 
limit k <^ uj <ti T as remarked above. 

Anticipating that the next-to-leading order correction 
will not dramatically change the leading order HTL dis- 
persion relations lu± (fc) which have a gap of order eT, we 
will first consider the case for which k ~ eT. As is clear 
from the above results, for generic k,Lu ^ eT the next-to- 
leading order correction to ReS(a;,fc) is of order T and 
gauge independent. Upon substituting ReS at this order 
into (20), we find that for the collective excitations of mo- 
menta k ~ eT the one-loop dispersion relations at next- 
to-leading order in T [i.e., 0{T)] are manifestly gauge in- 
dependent. This is one of the novel contributions of this 
article. Numerical analysis shows that the dispersion re- 
lations at next-to-leading order have similar features but 
move slightly above the leading order HTL results in the 
same momentum region. The differences increase with 
increasing gauge coupling constant, in agreement with 
the full numerical result found in Ref. [17]. 

The next-to-leading order correction to ImS in the mo- 
mentum region fc ^ eT is again of order T but gauge 
dependent. Whereas at next-to- leading order ImS has 
support above the light cone which presumably suggests 
finite damping rates 7 ^ e^T (up to the logarithm of the 
gauge coupling constant) for the collective excitations at 
one-loop order, nevertheless the gauge-dependent contri- 
bution to ImS that first appears at this order makes such 
interpretations doubtful. This result is not unexpected, 
as it is well-known that in high temperature gauge the- 
ories a consistent, gauge-independent calculation of the 
quasiparticlc damping rates requires a resummation of 
perturbation theory such as, e.g., the Braaten-Pisarski 
HTL resummation program [1]. 

We next consider the case for which k ^ eT. Ex- 
panding the fcrmion self-energy obtained above in the 
high-temperature limit [see (13), (15)-(17)] in powers of 
k/uj and keeping terms leading in k/uj and up to next- 
to-leading order in T, we obtain (fc <C tj ^ eT) 



and 



ReE(")(w,fc) 
ReE(^'(w,fc) 

ImS(°)(w,fc) 
ImS(i)(w,fc) 



e'T^ e'uj , T 

C — ~ 111 ■ 



8w 



•Stt^ 



k 

LU 



e^T^ e^uj T 



(26) 



e^T 
k 

LU 



(C-1) 



Se^T 
IGtt ' 



e^T ,, e^r 

DTT 127r 



(27) 



In the above expressions we have combined the gauge- 
independent and -dependent contributions. The next- 
to-leading order corrections to RcS and IniE are of 



order InT and T, respectively, and manifestly gauge- 
dependent. Upon substituting the above results into 
(20), we find that for the collective excitations of mo- 
menta k <C eT the one-loop dispersion relations and 
damping rates at next-to-leading order in T are gauge 
dependent. 

The effective fermion mass m as well as the damp- 
ing rate of collective excitations at rest 7(0) at next- 
to-leading order can be extracted by taking the long- 
wavelength limit k —f 0. The effective mass is determined 
by the following equation {m > 0) 



e^T^ e^m\ T 



+ f 
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(28) 



which agrees with the result found in the real-time for- 
malism [18]. Eq. (28) indicates that the one- loop effective 
mass at next-to-leading order in T is gauge dependent 
and that the gauge-dependent correction to m is of order 
e^T^ ln(l/e). The damping rate of collective excitations 
at rest is of order e^T but with a gauge-dependent con- 
tribution of the same order. 

Before ending this section we discuss briefly possible 
cancellation of the leading gauge dependence of the quasi- 
particlc poles at two-loop order. For simplicity we will 
consider the fc = case in QED. At two-loop order 
there are three one-particle- irreducible (IPI) diagrams 
contribute to the fermion self-energy, corresponding to 
corrections of the vertex, fermion self-energy and gauge 
boson polarization. However, because of the Ward iden- 
tity for the gauge boson polarization only the first two 
will give rise to gauge-dependent contributions. For the 
purpose of our discussion here, we only need the two- 
loop gauge-dependent contributions linear in (^—1). Pre- 
sumably, the gauge-dependent contributions quadratic in 
(^ — 1) will be subleading in T. 

To extract the leading gauge-dependent contribution 
at two-loop order, we focus on diagrams at the level of 
the propagator rather than the IPX self-energy diagrams. 
Specifically, we consider the one-particle-reduciblc (IPR) 
two-loop self-energy diagram. Using Y,^^' given in (26) 
and (27), we find for the real part (i) the leading gauge- 
independent contribution is of order e^T^/co^ correspond- 
ing to what will be resummed by the HTLs at two-loop 
order, and (ii) the leading gauge-dependent contribu- 
tion is of order (e^T^/w) ln(T/a;) arising from taking the 
gauge-dependent contribution in one of the two loops. 
Similarly, for the imaginary part we find both the lead- 
ing gauge-independent and -dependent contributions are 
of order e^T'^/uj^. 

Assuming that the leading gauge-dependent contribu- 
tions from the IPI diagrams have the same behavior as 
those from the IPR one, we obtain in fc ^ limit 



e2r2 .2_2 



7(0) = 
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,(29) 



where a(^) and 6(^) are ^-dependent coefficients yet to 
be determined by cxpHcit calculations. To solve for m 
and 7(0) at two-loop order, one can substitute into terms 
on the right-hand side in (29) the leading order result 
m ~ 0{eT). If o(^) and b{^) are such that the gauge- 
dependent contributions in (29) cancel, one finds the fol- 
lowing gauge-independent results 



e'T' 



O ( e^ In - 



7(0) ~ ©(e^T). (30) 



The above discussion seems to suggest that cancellation 
of the leading gauge dependence in m and 7(0) at two- 
loop order and next-to-leading odder in T is plausible. 
This situation is reminiscent of the recent proof that the 
truncated on-shell 2PI effective action has a controlled 
gauge dependence, with the explicit gauge-dependent 
terms always appearing at higher order [11-13]. 



We have discussed possible cancellation of gauge de- 
pendence at two-loop order in the case of the effective 
fermion mass and the damping rate for collective excita- 
tions at rest in QED. Our analysis suggests that cancel- 
lation of the leading gauge dependence at next-to-leading 
odder in T is plausible. A detailed study at two-loop or- 
der that allows us to verify such cancellation in hot gauge 
theories is the subject of further investigations. 
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IV. CONCLUSIONS 

In this article we have studied the gauge dependence of 
the fermion quasiparticle poles in hot gauge theories at 
one-loop order and next-to-leading order in T. We focus 
on the quasiparticle poles corresponding to soft collec- 
tive excitations of momenta k < eT <C T, with a view 
towards going beyond the leading order HTL results and 
resummations thereof. 

We have calculated the one-loop fermion self-energy 
(both the real and imaginary parts) in general covari- 
ant gauges up to next-to-leading order in the high- 
temperature expansion for which k^LO <^T. We find that 
whereas the next-to-leading order contributions to the 
imaginary part behaves like T, the behavior of those to 
the real part depends on the range of u as well as on their 
gauge dependence. The next-to-leading order gauge- 
independent contribution behaves like T for k,oj ^ eT, 
but becomes InT for fc <C w ~ eT. Nevertheless, the 
corresponding gauge-dependent contribution always be- 
haves like InT in the high-temperature limit. This anal- 
ysis allows us to study in detail the gauge dependence 
of the complex quasiparticle poles corresponding to soft 
collective excitations. 

For collective excitations of momenta k ^ eT we find 
that the dispersion relations at next-to-leading order in 
T are gauge independent. Numerical results shows that 
these next-to-leading order dispersion relations have sim- 
ilar features but move slightly above the leading order 
HTL ones in the same momentum region. However, the 
corresponding damping rates are gauge dependent, thus 
rendering the complex quasiparticle poles gauge depen- 
dent. For k <C eT and in the long- wavelength limit 
A: ^ 0, both the dispersion relations and the damping 
rates are gauge dependent. The gauge dependence of the 
position of the quasiparticle poles at one-loop order and 
next-to-leading order in T signals the need for resumma- 
tions of perturbation theory. 



APPENDIX A: VACUUM CONTRIBUTION TO 
THE FERMION SELF-ENERGY 

In this appendix we calculate the vacuum contribution 
to the one-loop fermion self-energy for the sake of com- 
pleteness. As a by-product, we show that the position 
of the singularities of the fermion propagator in vacuum 
are gauge independent at one-loop order. 

The vacuum contribution to the imaginary part of the 
self-energy ImE is extracted from (12) and (14) by ne- 
glecting the thermal distribution functions. The angular 
integration can be done analytically in the same manner 
as that in the finite-temperature contribution, and the 
remaining radial integration over momentum is elemen- 
tary. We find 



ImS™(S^±(a>,fc) = sgn(^) — (c. T k)9{L0^ ~ k^), 
lmJ:i%{cj, fc) = (e - 1) sgn(c.)-^(c. T k)9{cj' - k^), 
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where S± are defined in (21). We note that the gauge- 
independent and -dependent contribution can be com- 
bined together to yield a single gauge-dependent term 
proportional to ^. As expected, IniE has support only 
above the light cone corresponding to the usual two- 
particle cuts. 

The real part obtained from the imaginary one through 
the dispersive representation (8), is ultraviolet divergent 
and hence has to be regularized by an ultraviolet fre- 
quency cutoff. Introducing counter terms with the renor- 
malization condition ReS™^(/i) = at some arbitrary 
renormalization point /i, we obtain 



ReE^^"(w,fc) 



if! 
167r- 
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where the gauge-independent and -dependent contribu- 
tions have been combine together. 

Upon substituting the above results into (20), we find 
that the putative real poles of the fermion propagator in 
vacuum is determined by 



[u^k) 



1 
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(A3) 



Because of the infrared logarithmic divergence at thresh- 
old associated with the emission of soft gauge boson, 
w ^ ±fc are no longer isolated poles but the endpoints of 
logarithmic branch cuts. Nevertheless, we clear see that 
the position of the singularities (in this case logarithmic 



branch points and branch cuts) of the fermion propagator 
in vacuum are gauge independent at one-loop order. 



If one adds the vacuum and the finite-temperature con- 
tributions, one finds that the real part of the vacuum 
contribution can be exactly combined with the In T term 
in the real part of the finite-temperature contribution to 
yield a single term proportional to ln(T//i). This gen- 
eralizes the fc = case found in Ref. [18] to the case of 
fc,a; <C T and is the fermionic counterpart of the results 
found in Rcfs. [4, 20, 21]. However, no similarly combi- 
nation can be found in the corresponding imaginary part 
of the fermion self-energy. 
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